Entropic particle transport: higher order corrections to the Fick- Jacobs diffusion 
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Transport of point-size Brownian particles under the influence of a constant and uniform force field 
through a three-dimensional channel with smoothly varying periodic cross-section is investigated. 
Here, we employ an asymptotic analysis in the ratio between the difference of the widest and the most 
narrow constriction divided through the period length of the channel geometry. We demonstrate that 
the leading order term is equivalent to the Fick-Jacobs approximation. By use of the higher order 
corrections to the probability density we derive an expression for the spatially dependent diffusion 
coefficient D(x) which substitutes the constant diffusion coefficient present in the common Fick- 
Jacobs equation. In addition, we show that in the diffusion dominated regime the average transport 
velocity is obtained as the product of the zeroth-order Fick-Jacobs result and the expectation value 
of the spatially dependent diffusion coefficient (D(x)). The analytic findings are corroborated with 
the precise numerical results of a finite element calculation of the Smoluchowski diffusive particle 
dynamics occurring in a reflection symmetric sinusoidal-shaped channel. 

PACS numbers: 05.60. Cd, 05.40.Jc, 02.50.Ey, 51.20.+d 



I. INTRODUCTION 

The transport of large molecules and small particles 
that are geometrically confined within pores, channels or 
other quasi-one-dimensional systems attracted attention 
in the last decade. This activity stems from the prof- 
itableness for shape and size selective catalysis, particle 
separation and the dynamical characterization of poly- 
mers during their translocation I n particular, the 
latter theme which aims at the experimental determi- 
nation of the structural properties and the amino acid 
sequence in DNA or RNA when they pass through nar- 
row openings or the so-called bottlenecks, comprises chal- 
lenges for technical developments of nanoscaled channel 
structures 

Along with the progress of the experimental techniques 
the problem of particle transport through corrugated 
channel structures containing narrow openings and bot- 
tlenecks has give rise to recent theoretical activities to 
study diffusion dynamics occurring in such geometries 
Previous studies by Jacobs Q and Zwanzig [To| 
ignited a revival of doing research in this topic. The 
so-called Fick-Jacobs approach [To|. accounts for the 
elimination of transverse stochastic degrees of freedom 
by assuming a fast equilibration in those transverse di- 
rections @. The theme found its application for particle 
transport through periodic channel structures and 
designed single nanopores fl2j exhibiting smoothly vary- 
ing side walls. Several aspects of driven motion in pres- 
ence of applied external force fields and the quality of 
the Fick-Jacobs approach in presence of an applied force 
field in corrugated structures has been the focus of recent 
studies [H-Il|. 
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Beyond the Fick-Jacobs (FJ) approach, which is suit- 
ably applied to channel geometries with smoothly vary- 
ing side walls, there exist yet other methods for de- 
scribing the transport through v arying c hannel structures 
like cylindrical septate channels [19l - [2l| , tubes formed by 
spherical compartments [22l [23j or channels containing 
abrupt changes of cross diameters [24 1 . 

Our objective with this work is to provide a system- 
atic treatment by using a series expansion in terms of a 
smallness parameter which specifies the channel corru- 
gation for biased particle transport proceeding along an 
extended, three-dimensional periodic, reflection symmet- 
ric channel for which the original, commonly employed 
(lowest order) Fick-Jacobs approach fails because of ex- 
treme bending of the channel's side walls. 

In Sec.|ll]we introduce the model system: a Brownian 
particle in a confined channel geometry with reflection 
symmetric, irregular boundaries. The central findings, 
namely the analytic expressions for the probability den- 
sity and the average transport characteristics are pre- 
sented in Sec. IIHI In Sec. IIVI we employ our analytical 
results to a specific channel configuration consisting of si- 
nusoidally varying side walls. Section fVl summarizes our 
findings. 



II. TRANSPORT IN CONFINED STRUCTURES 

Generic mass transport through confined structures 
such as irregular pores and channels occurs due to the 
combination of molecular diffusion, as quantified by the 
molecular diffusivity D, and passive transport arising ei- 
ther from different particle concentrations maintained 
at the ends of the channel, an applied hydrodynamic 
velocity field or an external, force generating potential 
U(x, y, z). Here, we concentrate on constant force-driven 
transport where particles of dilute concentration (i.e. in- 



2 



A 




Figure 1: Sketch of a segment of a reflection-symmetric sinu- 
soidally varying channel that is confining the motion of the 
overdamped, point-like Brownian particle. The periodicity of 
the channel structures is L, the height AH , the minimal and 
maximal channel widths are Aw and Af2, respectively. The 
size of an unit-cell is indicated with the dashed lines. 



teraction effects can safely be neglected) are subjected to 
a fixed external force with magnitude F acting along the 
longitudinal direction of the channel e x , i.e., U(x, y, z) — 
—Fx. The overdamped single Brownian particle then 
budges in a three-dimensional periodic channel geometry 
of period L, constant height AH , and periodically vary- 
ing transverse width. A sketch of a segment of the chan- 
nel is depicted in Fig. [TJ The shape of the side walls are 
described by the two boundary functions uj±(x). As we 
restrict ourselves to reflection-symmetric confinements in 
y-direction, we set u>±(x) = ±lo(x). 

The evolution of the probability density P (q, t) of 
finding the particle at the local position q = (x, y, z) T 
at time t is governed by the three-dimensional Smolu- 
chowski equation [H|,[26|, i.e., 



&P(q,i)+V q . J(q,t) = 0. 



(la) 



where 



J(q,t) = -P(q,t)e x -^V q P(q,i) (lb) 

is the probability current of the probability density 
P (q, t) . The force strength acting on the Brownian 
particle is denoted by F, r\ is the friction coefficient, 
while the Boltzmann constant is ks and T refers to 
the environmental temperature. Because of the impen- 
etrability of the channel walls the probability current 
J (q, t) = (J x , J y ', J Z ) T is subjected to the no-flux bound- 
ary condition, reading 



J (q, t) ■ n = , Vq € channel wall . 



(2) 



n denotes the out-pointing normal vector at the channel 
walls. The probability density satisfies the normalization 
condition / - t _ n P(q, t) d 3 q = 1. as well as the period- 
icity condition P(x + mL, y, z,t) — P(x, y, z,t) , Vm 6 Z. 
In the long time limit the stationary probability den- 
sity is defined as P s t (q) := limt_ ) . 00 P (q, t) . Analo- 
gously, the stationary probability current reads J s t (q) := 
limt^oo J (q,t). 

The key quantities of particle transport through such 
periodic channels are the average particle velocity (q) 



and the effective diffusivity D e g. The latter is given by 

(x*(t))-{x(t)f 



-D e ff = lim 



2t 



(3) 



and can be calculated by means of the stationary proba- 
bility density P s t (q) using an established method taken 
from Ref. [27|. Once P s t(q) is known, the mean particle 
velocity of Brownian particles can be computed by 



lim 



<q(*)> 



J st (q) d 3 q. 



(4) 



We next introduce dimensionless variables. In doing 
so, we measure longitudinal length and height as x = x/L 
and z = z/L, respectively. For the rescaling of the y- 
coordinate, we introduce the dimensionless aspect pa- 
rameter s: This is the difference of the widest cross- 
section of the channel, i.e. AO, and the most narrow 
constriction at the bottleneck, i.e. Aw, in units of the 
period length, yielding 



£ = 



(An - Alu) 
L 



(5) 



The dimensionless value of e characterizes the deviation 
of the boundary from the straight channel which amounts 
to e — 0. Following the reasoning in Ref. [28j], we next 
measure, for the case of finite corrugation e / 0, the 
transverse length y in units of eL, i.e. y = sLy and, like- 
wise, the boundary functions h±(x) = uj±(x)/(eL). Time 
is measured in units of r = L 2 rj/ '(/eg T) which is twice 
the time the particle assumes to overcome diffusively, at 
zero bias F = 0, the distance L, i.e. t = t/r. The poten- 
tial energy is rescaled by the thermal energy k^T, i.e., for 
the considered situation with a constant force component 
in channel direction: U = —Fx/jknT) = — fx, with the 
dimensionless force magnitude [lj], LL4[ : 



/ = 



FL 
k^T 



(6) 



The dimensionless forcing parameter / is given as the ra- 
tio of the work F L done on the particle when dragged by 
the constant force F along a distance of the period length 
L divided by the thermal energy ksT. Note, that for an 
adjustment of a certain value of / in an experimental 
setup one can modify either the force strength F or the 
temperature T. After scaling the probability distribution 
reads P (q, t) — e L 3 P (q, t) , respectively, the probabil- 
ity current is given by J (q,<) — t L 2 (eJ x , J v ,eJ z ) T . In 
the following, we shall omit the overbar in our notation. 

In dimensionless units, the Smoluchowski equation, cf. 
Eqs. (HJ, reads: 



ftP(q,i)+V q - J(q,i) = 0, 



(7a) 



where V q = (d x , ~d v , d z ) and 

J(q,t)= /P(q,*)e x -V q P(q,t) . (7b) 
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At steady state, Eq. (f7a|) becomes: 

s 2 d x j: t + d y jy t + e 2 d z j: t = o. 



(8) 



Because (i) the dynamics in z-direction is decoupled from 
the dynamics in x and y-direction and (ii) the shape of 
the lower and upper boundary depends neither on x nor 
on y, the separation ansatz P st (x,y,z) — p st (x,y) £(z) 
and the boundary condition 

J* t =0 , at z = and z = AH/L , (9) 

results in a non-trivial solution for £(z) for J s z t (q) = 
everywhere within the channel. For the considered sit- 
uation, i.e. there is only a constant force acting in re- 
direction, the form function £(z) equals the inverse of 
the dimensionless channel height, i.e. C = L/ AH. Note, 
that the presented separation technique can also be ap- 
plied for more complex forcing scenarios. Assuming a 
general potential landscape U(x,y,z) = V(x,y) + W(z) 
defined within the channel, the used separation ansatz 
for the stationary solution results in 

e -W{z) 

P st (x,y,z) = Pst (x,y) ■ • ( 10 ) 

o 

Consequently, this allows a reduction of the problem's 
dimensionality from 3D to 2D: 



(11) 



Note, that the 2D transport problem was investigated 
in symmetric [ll|, l^ - fla l28Tj30j and asymmetric [H, 
|24| channels. For an arbitrary dimensionless chan- 
nel geometry h± (x) the outwards pointing normal vec- 
tor at the perpendicular side walls is given by n = 

(jfh±(x), ±1, 0^ / ^Jl + h'±(x) 2 with the prime denot- 
ing the differentiation with respect to x. Therefore, the 
no-flux boundary condition Eq. ([2]) can be written as 



e 2 h' ± (x)J* t =J* t , Vyeh ± (x). 



(12) 



Note that even in the case of a more general substrate 
potential given by U(q) = V(x, y)+W(z) the 2D problem 
Eq. (fTTj) does not dependent on the potential W{z). 

Finally, we define the marginal one-dimensional prob- 
ability density in force direction p st (x) as follows 



h + (x) AH/L 



Pst (x) 



dy 



dzP st (x,y,z) 



(13) 



h_(x) 



III. ASYMPTOTIC ANALYSIS 

We apply the asymptotic analysis [HI, Hi], [HJ to the 
problem stated by Eq. (fTTj) and Eq. (|T2"j) . In doing so, we 



use for the stationary probability density p st {x,y) (the 
index st will be omitted in the following) the ansatz 



P(x,y) = ^2e 2n p n (x,y) , 



71 = 

and for the probability flux 



J(x,y) = Y,£ 2nj n(x,y) 



(14) 



(15) 



n=0 



in the form of a formal perturbation series in even orders 
of the parameter e. Substituting these expressions into 
Eq. (dU, we find 



= d y jy(x, y) + Y, e 2 " {8*JZ-i(x, V) + dyJ^(x, y)} , 

(16a) 



n=l 



and the no-flux boundary condition at the channel walls 
Eq. (JTJJ) turns into 

OO 

0= -4{x,y) + Y i ^ n {h ± {x)JZ-i{x,y)-Jl{x,y)) . 

n=l 

(16b) 

Each order p n has to obey the periodic boundary condi- 
tion p n (x + m, y) = p n (x, y) , Vm e Z and p(x, y) has to 
be normalized for every value of e. 

Consequently, the average particle velocity is given by 

OO 

(*) =(x) Q + J2e 2n {/ (Pn(x,y)) XtV - (d xPn (x,y)) x y } . 

71=1 

(17) 

In Eq. (|17l) . the average of an arbitrary function k(x,y) 
is defined as the integral over the cross-section in y 
and over one period divided by the period length which 
is one in the considered scaling, i.e. (k(x, y)) x y — 

Jq dx J_M x )dy k(x,y). In Sec. IIII Al we demonstrate 
that the zeroth order of the perturbation series expansion 
coincides with the Fick- Jacobs equation [H,[l(|. Referring 
to [HI HI] an expression for the average velocity (x) is 
known. Moreover, in Sec. IIIIB1 the higher orders of the 
probability density are derived. Using those results we 
are able to obtain corrections, see in Sec. IIII C[ to the 
average velocity beyond the zeroth order Fick- Jacobs ap- 
proximation presented in the next section. 

A. Zeroth Order: the Fick-Jacobs equation 

For the zeroth order, Eqs. (fl"6|) read 



d y J y Q {x,y) = -d y e- v ^d y (e v ^p (x,y)) =0, 

(18a) 
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supplemented with the corresponding no-flux boundary 
condition 



Consequently, 



4(x,y) =0, Vy ewaU. 



P0 (x,y)= g{x)e- y ^ 



(18b) 



(19) 



where g(x) is an unknown function which has to be de- 
termined from the second order O (e 2 ) balance given by 
Eq. (|L6a|) . Integrating the latter over the cross-section in 
y and taking the no-flux boundary conditions Eq. (|16b[) 
into account, one obtains 



0=d x 



-A(x) 



9'{x) 



(20) 



where the effective potential A(x) is explicitly given by 

(21) 



+h{x) 

-A(x) = f dye ~V( X ,y) 
-h(x) 



For the problem at hand, i.e. for V{x, y) — — fx, as well 
for potentials where x enters only linearly and where x 
is not multiplicatively coupled to the other spatial co- 
ordinates [lo, H3, HH the stationary probability density 
within the zeroth order reads 



e -V(x,y) X j e A(x') dx/ 

X 

~T x+i 

/ dxe-^ J e A ( x ')dx' 

x 



p (x,y) =e- v ^g{x) 



In addition, the marginal probability density Eq. (113 
becomes 



(22) 



Po(x) 



-A(x) 



(23) 



Expressing next g(x) by po(x), see Eq. (|2"0"|) . then yields 
the celebrated stationary Fick- Jacobs equation 



= d x (e- A ^d x e A ^p o (x) 



(24) 



derived previously in Ref. || [T(| H||- Thus, we find 
the result that the leading order term of the asymptotic 
analysis is equivalent to the FJ-equation. Please note 
that the differential equation determining the unknown 
function g(x), cf. Eq. ([20)) . is the same for the dynam- 
ics of a Brownian particle evolving in an energetic po- 
tential V en (x,y) leading to a confinement in y-direction, 
with the natural boundary conditions J%(x,y = ±oo) = 
[l(J HH- Therefore, in zeroth order and for the 
given scaling, an appropriately chosen confining energetic 
potential V cn (x,y) obeying J^dy exp(-V cn (x,y)) = 

J^),dyexp(—V(x,y)) results in the same transport 



characteristics as those induced by the confining chan- 
nel with the boundary functions h± (x) [37j • 

The average particle current is calculated by i nteg rat- 
ing the probability flux Jg over the unit-cell [33] , 1381 ] 

1 h + (x) 

( i (/)>o = dx / dyj*(x,y) 



(25) 



ft_(x) 

1 - e~f 



1 x + l 

Jdxe'^ J e A (*')dx' 

x 



In the spirit of linear response theory, the mobility in 
units of the free mobility 1/r] is defined by the ratio of 
the mean particle current Eq. (|25[) and the applied force 
/ yielding 



Wo (/) 



<*(/)>„ 



(26) 



B. Higher order contributions to the Fick- Jacob 
equation 

We next address the higher order corrections p n (x,y) 
of the probability density. According to Eq. (|16a.|) . one 
needs to iteratively solve 



d%p n {x,y) =£pn-i(x,y) , n>i, 



(27) 



under consideration of the boundary condition Eq. (|16bl) . 
In Eq. |27p , we make use of the operator £, reading £ = 
{f d x — d x ). Applied n-times yields the expression 



fc=0 v 7 



(28) 



Each solution of the second order partial differential 
equation Eq. (|2T)1 possesses two integration constants 
d n ^i and c?„.2- The first one, d n> \, is determined by the 
no-flux boundary condition Eq. (I16b[) while the second 
provides the normalization condition {p(x,y)) x = 1. In 
what follows, we use the normalization constant of the 
probability density p{x, y) via the zeroth order (po(x, y)} . 
As a consequence, we have 

1 h + (x) 

(Po(x,y)) xy = J dx J dyp (x,y)=l, (29a) 

h-(x) 

(jPn(x,y)) Xjy =0, Vn>l, (29b) 
with the constraint that 

h+(x) 

dyp n (x,y)^0, Vn>l, (29c) 



h_(x) 
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in order to prevent that the marginal probability density 
Eq. (fig)) equals the FJ results, cf. Eq. (j2"3"]) . for an arbi- 
trary value of e, i.e. p(x) — po(x). Further, we have to 
emphasize that the centered functions 



Pn(x,y) 



Pn (X,V) ~ (Pn (x,y)) 



for n > 1 , (30) 



are no probability densities anymore because they can 
assume negative values for a given x and y. The calcu- 
lation of the average particle velocity Eq. (flTf simplifies 
to 



(x) = (x) n - 



DO 

n=0 



.2;i 



{d x p n (x,y)) xy 



(31) 



We find that the average particle current is composed of 
(i) the Fick- Jacob result (x) , cf. Eq. (T2~5l) , and (ii) be- 
comes corrected by the sum of the averaged derivatives of 
the higher orders p n (x, y). One immediately notices that 
the second integration constant d n ,2 does not influence 
the result for the average particle velocity Eq. (13T1) . 
For the first order correction, the determining equation 



dyPi(x,y) = £p (x,y) = —^-d x 



h{x) 



(32) 



and after integrating twice over y, we obtain 



with the integration constants for the n-th order 




d n , 2 = - / dx ^£"- fe d M 



k=l 



dx£ n p {x,y) 



h 2n+l 

(2n+T)! 



^2(n-fc)+2 

(2 (n-k) + 2)\ 



dx h{x) . 



(35b) 



As expected, for a reflection symmetric channel in y- 
direction each order p n (x,y) results as well in a reflec- 
tion symmetric function. The latter consists of a term 
proportional to even powers in y and in addition of a 
sum of odd powers of \y\, caused by the no-flux bound- 
ary conditions. Since each integration constant d nt i(x) 
with n > 1 is determined by the probability current of 
the previous order, every order p n (x,y) is proportional 
to the average current of the zeroth order (x) Q . Con- 
sequently, the 2D probability density equals the zeroth 
order p{x, y) = po(x,y) = const for all values of e in 
absence of an external force / = 0. Further, it follows 
that the average particle current Eq. (T3"Tj) scales with the 
average particle current obtained from the Fick- Jacobs 



formalism 



for all values of e. 



h'(x) 
h 2 (x) 



y_ 

2! 



(33) 



C. Spatially diffusion coefficient D(a 



Hereby, as previously requested above, the first integra- 
tion constant di t i(x) is set to in order to fulfill the 
no-flux boundary condition, and the second must pro- 
vide the normalization condition Eq. (|29b[) . i.e. d\$ = 0. 
Consequently, the first correction to the probability den- 
sity becomes positive if the confinement is constricting, 
i.e. for h'(x) < and (x) ^ 0. In contrast, the prob- 
ability density becomes less in unbolting regions of the 
confinement, i.e. for h'(x) > 0. Please note, that the first 
order correction scales linearly with the average particle 
current (x) Q . Overall, the break of spatial symmetry ob- 
served within numerical simulations in previous works 
El 39] is reproduced by this very first order correction. 
Particularly, with increasing forcing, the probability for 
finding a particle close to the constricting part of the 

confinement increases, cf. Ref. mm. 

Upon recursively solving, we obtain for the higher or- 
der corrections n > 1 as 



2)1 



y 

Pn (X, y) = £ n P(){x, y)^J+ d n,2- 



J2& n - k d kll ( X ) 



|2(n-fe) + l 



(34) 



k=l 



(2 (n-fc) + l)! 



With Sec. IIII Al we could show that the dynamics 
of Brownian particles in confined structures can be de- 
scribed approximatively by the FJ-equation, cf. Eq. (j2"4")l . 
Zwanzig [Hj obtained this ID equation from the full 2D 
Smoluchowski equation upon eliminating the transverse 
degree of freedom. This approximation neglects the in- 
fluence of relaxation dynamics in transverse direction, 
supposing that it is infinitely fast. In a more detailed 
view, we have to notice that diffusing particles pile up, 
or miss, at the curved wall if the channel is getting nar- 
rower or wider as they can flow out from/ or towards the 
wall in y direction only at finite time. These effects are 
described by the higher expansion orders p n (x,y) pre- 
sented in Eq. (|34[) . In the following, we aim at deriving 
a dynamical equation of Smoluchowski- type, but with 
a diffusion coefficient that depends on the longitudinal 
channel coordinate x. 

The concept of a spatially dependent diffusion coef- 
ficient D(x) was introduced by Zwanzig [10( and sub- 
sequently supported by the study of Reguera and Rubi 
[361 ] . The main idea is to combine all marginal higher 
orders corrections p n (x) into a one-dimensional func- 
tion D(x), effectively acting on the marginal probabil- 
ity density p(x). In Ref. [361 ] the corrected stationary 
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FJ-equation has the form 







-d x J x (x) 



d_ 

d x 



DMe^W^WpW , (36) 
ox 



and was derived therein within the framework of meso- 
scopic non-equilibrium thermodynamics. Kalinay and 
Percus [HI used a rigorous mapping of the 2D diffusion 
equation onto the reduced dimension and derived an ex- 
pansion of the diffusion coefficient D(x), which represents 
corrections to the FJ-equation. 

In this spirit we now determine the spatial dependent 
diffusion coefficient D(x) based on the presented results 
for the perturbation series expansion Eq. We con- 

centrate on the limit of small force strengths |/| <C 1, so 
that diffusion is the dominating process. Integrating the 
2D stationary Smoluchowski equation Eq. (jlll) over the 
cross-section in y, and respecting the no-flux boundary 
condition Eq. (|12p . one derives an alternative definition 
of the marginal probability current J x (x), equivalent to 
Eq. (EH): 



-J x (x) =D(x)h(x)d x 



h(x) 



h(x) 



dxp{x, y)dy . 



-h(x) 



(37) 



The second equality determines the sought-after spatial 
dependent diffusion coefficient D(x). Note, that D(x) 
is solely determined by derivatives of p(x,y) and p{x). 
Hence, it plays no role whether one uses the original ex- 
pansion terms defined by Eq. (|I6a[) or the centered ones, 
given by Eq. (|30[) . In compliance with Ref. we make 
the ansatz that all but the first derivative of the bound- 
ary function h(x) are negligible. Then, the integration 
constants d ni i(x) equal as they can been shown to be 
proportional to higher derivatives of h(x). Moreover, in 
the limit |/| <C I, the n-times applied operator £, cf. 



Eq. <J2SJ), simplifies to £" = (-1)* 



£> = (i) (-I)"(2n-l)! 



/\2n-l 



2h 2r 



Moreover 



0{h"(x)). (38) 



Inserting the probability densities into Eq. (l37j) , one finds 
that 



D(x)=J2e 2n (-l) 



(h 



l\2n 



n=0 

arctan (eh'(x)) 
' eh'(x) 



2n + 1 



0(h"(x)) 



(39) 



for the spatially dependent diffusion coefficient D(x) in 
the diffusion dominated regime, i.e. when |/| <C I . Note, 
that this expression for D(x) was obtained previously by 
Kalinay and Percus (l5j within a quite different expan- 
sion approach. 

In what follows, we evaluate the average particle cur- 
rent Eq. (|3ip by means of the spatially dependent dif- 
fusion coefficient D(x). According to lim/_ i .o (x) = 



dxJ x (x), it follows that in the small force limit the 
mean particle velocity is proportional to the expecta- 
tion value of the spatially dependent diffusion coefficient, 
yielding the main finding 



fm {*(/)) = Jim <±(/)>o (D{x)) x + O (h"(x)) 

,. wj. N v / arctan (eh' (x)) 
^?o W) H eh>{x) 



(40) 



In Eq. (|40|) . the average of an arbitrary function k(x) 
is defined as the integral over one period divided by the 
period length which is one in the considered scaling, i.e., 
(k(x)) x = J Q k(x) dx. In the linear response limit, i.e. 
for | /| <C 1, the Sutherland- Einstein relation emerges 
Ref. 0|4l[, reading in dimensionless units: 



limM/) = lim £>eff(/), 



(41) 



the effective diffusion coefficient D e g is determined by the 
mobility fi = lim/^o ( x (f)} If- Consequently, if the av- 
erage current (i(/)) (or the effective diffusion coefficient 
D® s (f)) are known in the zeroth order, the higher order 
corrections to both quantities can be obtained according 
to Eq. flM]). 



IV. APPLICATION OF THE THEORY TO A 
SINUSOIDALLY SHAPED CHANNEL 

In the following we validate the obtained analytic pre- 
dictions Eq. (|40[) with precise numerical simulations con- 
cerning one single point-like Brownian particle mo ving 
with a corrugated sinusoidally-shaped geometry Ref. [13l . 
[T3 |. The dimensionless boundary function h(x) reads 

h±{x) = ±h(x) = ±i (^^+sm(2nx) \ , (42) 

and is illustrated in Fig. [TJ Please note, that in absence 
of the scaling each channel geometry is determined by the 
period L, the maximum width Ai7, and the width at the 
bottleneck Atu. Upon scaling all lengths are measured 
in units of the period L. Consequently, the parameter 
5£l denotes the ratio of the maximum width Af2 and 
the period L, viz., Sfl = Afl/L. Equivalcntly, it holds 
that 5uj = Auj/L. Within this scaling, the period of the 
channel equals one. 

In addition, one notices that the dimensionless bound- 
ary function h(x) is solely governed by the aspect ratio 
of the minimal and maximal channel width S = Sfl/5oj. 
Obviously different realizations of channel geometries can 
possess the same value of 5. The number of orders have 
to taken into account in the perturbation series Eq. (fl"4"| . 
respectively, the applicability of the Fick- Jacob approach 
to the problem, depends only on the value of the slope 
parameter e — Sfl (1 — S) for a given aspect ratio S. For 
clarity, the impact of the maximum S£l and minimum 
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Figure 2: (Color online) Schematic sketch of the dependence 
of the expansion parameter e = S£l — Sui and the aspect ra- 
tio 5 = 5u/8S} on the maximum width 5Q, respectively, the 
width at the bottleneck 5lj in units of the period L. The 
dashed lines correspond to 5 — 1,0.5, 0.25 (from above) while 
the colored areas illustrate pairs of (SQ, 5uj) where e < 0.1 
(blue, circles), e < 1 (red, triangles), e < 5 (green, dots), and 
e > 5 (yellow,plus signs). 



width Suj on the expansion parameter e, respectively, the 
aspect ratio S is illustrated in Fig.O 

According to the Sutherland-Einstein relation Eq. (|4Tj) 
the mobility equals the effective diffusion coefficient (in 
the dimensionless units) for / <C 1 Consequently, it 
is sufficient to discuss the behavior of the mobility //(/). 
Referring to Sec. IIII CI the higher order corrections to 
the mobility are given by the product of the F J- result and 
the expectation value of the spatially dependent diffusion 
coefficient D(x), see Eq. (|40|) . 

First, we obtain the mobility /j,q within the zeroth- 
order (Fick- Jacobs approximation) . In the diffusion dom- 
inated regime, the analytic expression for the mobility 
within the FJ-approach, cf. Eqs. (|25|) and ([26]) . simpli- 
fies to the Lifson- Jackson formula fijl, 



:= lim /Lt (/) = 

(h(x)) 



1 



l 

h(x) 



lim D eS (f) . 
/->o 



(43) 



For the exemplarily considered channel geometry 
Eq. (|42|) the mobility attains the asymptotic value 



limM/) 



2VS _ 

1 + 5 ~ 2-e/Sn 



(44) 



One notices that in the diffusion dominated regime | / 1 <C 
1 the mobility of one single particle is determined only by 
the geometry - more precisely by the aspect ratio S. In 
the limit of vanishing bottleneck width, i.e. S — > 0, the 
mobility tends to 0. In contrast, for straight channels 
corresponding to 5 = 1, i.e. e = 0, the mobility equals 
its free value which is one in the considered scaling. 

Evaluating the period-averaged value of D(x), i.e., con- 
sidering all higher order corrections apart from than scal- 
ing with higher derivatives of the boundary function h(x), 



Figure 3: (Color online) Comparison of the analytic theory 
versus precise numerics (in dimensionless units): The mobil- 
ity and the effective diffusion constant for a Brownian particle 
moving inside a channel confinement are depicted as function 
of the ratio of slope parameter e and maximal channel width 
SQ for different values 8Q — 0.1,1,2,5 and bias / = 10~ 3 
(corresponding to the diffusion dominated regime). The sym- 
bols correspond to the numerical obtained mobility (triangles) 
and the effective diffusion coefficient (circles). The lines cor- 
respond to analytic higher order result, cf. Eq. (|45[) . The 
zeroth order - Fick- Jacobs results given by Eq. (I44|) collapse 
to a single curve hidden by the solid line. 



we obtain from Eq. (|4"0")l : 



M := ]™M(/) = Mo {D{x)) 



4^/1 - e/Sn asinh(7re/2) 
2-e/SQ Ve 



(45) 



for the mobility /i and the effective diffusion coefficient 
D e s in units of its free values, respectively. 

In Fig. [3J we depict the dependence of the /j,(f) (trian- 
gles) and D c g(f) (circles) on the slope parameter e for 
/ = 10~ 3 . The numerical results are obtained by solv- 
ing the stationary Smoluchowski equation Eq. (|lap using 
finite element method [43| and subsequently calculating 
the average particle current Eq. In order to deter- 
mined the effective diffusion coefficient D c g(f), one has 
to solve numerically the reaction-diffusion equation for 
the B-field [13, [H)]. Note, that the numerical results for 
the effective diffusion coefficient D c s(f ) and the mobil- 
ity /j,(f) coincide for all values of e, thus corroborating 
the Sutherland-Einstein relation. In addition, the Fick- 
Jacobs result, given by Eq. (|4"4"| . and the higher order 
result, see Eq. (|45[) . are depicted in Fig. [3] 

For the case of smoothly varying channel geometry, i.e. 
Sfl -C 1, the analytic expressions are in excellent agree- 
ment with the numerics, indicating the applicability of 
the Fick- Jacobs approach. As long as the extension of 
the bulges of the channel structures is small compared to 
the periodicity, sufficiently fast transversal equilibration, 
which serves as fundamental ingredient for the validity of 
the Fick- Jacobs approximation is taking place. In virtue 
of Eq. ([5]), the slope parameter is defined by e = Sfl — 5u> 



8 



and hence the maximal value of e equals 5£l, see Fig. [2] 
Consequently the influence of the higher expansion orders 
e 2n (d x p n (x, y)) on the average velocity Eq. (|3"l"j) and mo- 
bility, respectively, becomes negligible if the maximum 
channel's width 8fl is small. 

With increasing maximum width the difference be- 
tween the FJ-result and the numerics is growing. Specif- 
ically, the FJ-approximation resulting in Eq. (|44p over- 
estimates the mobility \i and the effective diffusion co- 
efficient D c ff- The higher order corrections need to be 
included and consequently provide a good agreement for 
a wide range of e-values for maximum widths SQ on the 
scale to the length of the channel, i.e. <JO ~ 1, see the 
dotted line in Fig. [3] Upon further increasing the max- 
imum width SD, diminishes the range of applicability of 
the derived higher order corrections. This is due to the 
neglect of the higher derivatives of the boundary func- 
tion h(x). Put differently, the higher derivatives of h(x) 
become significant for Sfl^> 1. 

V. SUMMARY AND CONCLUSION 

In summary, we have considered the transport of point- 
size Brownian particles under the influence of a constant 
and uniform force field through a three-dimensional chan- 
nel. The latter exhibits a constant height and periodi- 
cally varying side walls. 

We have presented a systematic treatment of parti- 
cle transport by using a series expansion of the station- 
ary probability density in terms of a smallness parameter 
which specifies the corrugation of the channel walls. In 
particular, it turns out that the leading order term of 
the series expansion is equivalent to the well-established 
Fick- Jacobs approach 0, The higher order correc- 
tions to the probability density become significant for 
extreme bending of the channel's side walls. Analytic re- 
sults for each order of the perturbation series have been 



derived. Interestingly, within the presented perturba- 
tion theory, all higher order corrections to the stationary 
probability distribution and the average particle current 
scale with the average particle current obtained from the 
Fick- Jacobs formalism. Moreover, by using the higher 
order corrections we have derived an expression for the 
spatially dependent diffusion coefficient D(x) which sub- 
stitutes the constant diffusion coefficient present in the 
common Fick- Jacobs equation. Accordingly, in the lin- 
ear response regime, i.e. for small forcing |/| <C 1, the 
mean particle velocity is then given by the product of the 
average particle current obtained from the Fick- Jacobs 
formalism and the expectation value of the spatially de- 
pendent diffusion coefficient D(x) . Moreover, due to the 
Sutherland-Einstein relation, the above statement also 
holds good for the effective diffusion coefficient. 

Finally, we have applied our analytic results to a spe- 
cific example, namely, the particle transport through a 
channel with sinusoidally varying side walls. We cor- 
roborate our theoretical predictions for the mobility and 
the effective diffusion coefficient with precise numerical 
results of a finite element calculation of the stationary 
Smoluchowski-equation. In conclusion, the considera- 
tion of the higher order corrections leads to a substantial 
improvement of the Fick- Jacobs-approach, which corre- 
sponds to the zeroth order in our perturbation analysis, 
towards more winding side walls of the channel. 
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